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Abstract: A feasibility study to extract the polarised gluon density from di-jet 
events at HERA in next-to-leading order is presented. It is shown that when tak- 
■ ing next-to-leading order effects into account the asymmetries at HERA remain 

measurable and sensitive to the polarised gluon distribution. They can be used to 
extract the polarised gluon density in the proton in the region 0.005 < x g < 0.4. 
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1 Introduction 

a: 



The origin of the spin in the proton is still a subject of much debate. Over the years it has 
been confirmed that the quarks, as measured in deep inelastic scattering, account for only 
30% of the proton spin. Next-to-leading order (NLO) QCD fits of structure function and semi- 
inclusive data suggest that the contribution of the gluon to the spin could be large |], |2|] . A first 
^ attempt of a direct measurement of the polarised gluon distribution, AG, using leading charged 
particles || is not in conflict with this suggestion. In general, it has been concluded that major 
progress in our understanding of the spin structure can be made with clear and unambiguous 
direct measurements of AG. Several experiments are planned to tackle this measurement [^|, |5|. 

It has been shown that the ep collider HERA, when both beams are polarised, could make 
an important contribution to the determination of AG(x), for a considerable x-range, where x is 
the momentum fraction of the proton carried by the gluon. A particularly sensitive method is to 
extract AG(x) from di-jet events in deep inelastic scattering. Feasibility studies of extracting 
the polarised gluon density AG(x) from di-jet events at HERA in leading order (LO) have 
been performed |], [7|, [J, the most detailed one was published in the proceedings of the 



workshop 'Future Physics with polarized beams at HERA' |J. The event generator PEPSI [10 
which includes hadronization of the final parton state, was used, followed by a simple detector 
simulation. PEPSI includes LO matrix elements for the QCD processes in the hadronic final 
state. A first estimate of higher order effects was obtained by including initial and final state 
unpolarised parton showers. 

Recently next-to-leading order polarised cross sections for di-jet production became available 
with the program MEPJET |]TT], [T!J. This allows to check NLO corrections to the LO cross 
section asymmetries. In this paper MEPJET will be used to find optimal event and di-jet 
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Figure 1: Feynman diagrams for the process of Boson-Gluon-Fusion (BGF) in LO (left) and a 
higher order process (right). 



selection cuts, which are a balance between statistical significance and analysing power of the 
asymmetries. 

Note that at LO the momentum fraction of the proton carried by the gluon, Xgi C3jI1 be 
directly calculated from the di-jet kinematics (see Fig. 1). We define 




calculated from the Bjorken-x, the four momentum transfer Q 2 and the invariant mass s of the 
di-jet system. At LO the variable Xj e t s is identical to x g , while at NLO Xj e t s is different from 
x g even at parton level e.g. due to gluon radiation processes, see Fig. |l|. In fact the relation 
x g > Xjets holds. 

Since the main purpose of this paper is to show the size of the asymmetries in NLO and 
the range in x g where one can expect information on AG, a simple method to relate the 
reconstructed Xjets variable from the jets to the true one will be used. The same technique to 
unfold the gluon from the data as done at LO || will be used. 

In the last chapter we show the potential of a future high energy ep collider, to extract AG 
in LO from di-jet events, using TESLA and HERA as an example. 



2 Di-jet selection using MEPJET 

We start with kinematic cuts close to the ones used in the published LO study || |9|: 

5 < Q 2 < 100 GeV 2 (1) 

y > 0.3 (2) 

Electron > 5 GeV (3) 

p{ ets > 5 GeV (4) 
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Table 1: Unpolarised and polarised di-jet cross-sections for different kinematic cuts, in LO and 
NLO. The cuts are the lower limits. The last two columns contain a measure for the analysing 
power (see text). 

Compared to @, [|, for the scattered electron and jet selection in this analysis, a pseudorapidity 
cut of \rj\ < 3.5 is used. An additional cut on the p J t ets > 5 GeV in the Breit frame is also 
applied. As before the cone jet algorithm was used with a cone size of 1 in azimuthal angle 
and pseudorapidity. For the calculations we used the structure function parametrisations of 
GRV |13[ for the unpolarised case, and Gehrmann-Stirling set A [[14]] for the polarised case. 



The cut on the square of the invariant mass of the two jets was varied in the range, s := 

Q 2 ( x jets/ X ~ I); as 

8 > 100,200,300,500 GeV 2 , 

where s is computed from the reconstructed jet quantities, and the variables are indicated in 
Fig. 1. Further cuts which were varied in the study are: 

p\ ets > 5, 7 GeV, 

Q 2 > 2,5,10,40 GeV 2 , 
y > 0.15,0.3, 

and finally the electron energy cut was lowered to E e = 3 GeV 2 . The results for LO and NLO, 
polarised and unpolarised cross section for exclusive di-jet production for several combinations 
of cuts, can be seen in Table 1. In the last column of this table the ratio between the expected 
overall asymmetry (Acr/a) and a quantity which is proportional to the expected statistical error 
(l/y/a) has been computed, and hence gives an idea of the sensitivity: a higher value means 
higher sensitivity. 
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One observes that for all tried scenarios the sensitivity decreases from LO to NLO. Further- 
more globally the corrections to the unpolarised cross sections become large for a high s cut, 
while the polarised cross section receives high corrections when a low s cut is used. For the final 
study we selected those scenarios where the correction to the polarised and unpolarised cross 
sections are less than 30% and a good compromise between the analysing power and statistics is 
found. This can be obtained by choosing s > 300 GeV 2 or s > 200 GeV 2 Ap t > 7 GeV, both of 
which have a large sensitivity to the gluon according to Table 1. The asymmetries as a function 
of Xj ets for these scenarios in LO and NLO are compared in Fig. 0. The values are reduced in 
NLO compared to LO, but still sufficiently large. The errors correspond to the statistical errors 
expected for 200 pb _1 (but assuming a 100% polarisation of the colliding beams). Table 1 also 
shows that lowering the electron energy requirement does not bring a significant improvement, 
despite the region of larger depolarisation factor included. We will therefore use the scenario 
of s > 200 GeV 2 A p t > 7 GeV and E e > 5 GeV in the following. 
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Figure 2: Comparison of the asymmetries versus Xj e t s in LO and NLO for the two selected cut 
scenarios (see text). 



Before moving towards extracting the gluon distribution we first compare di-jet asymmetries 
in LO predicted by MEPJET with PEPSI without parton showers, and find them to agree very 
well. This has been shown before ||, and has been confirmed here. In Table 2 the overall 
asymmetries are shown using different s cuts and event/jet selection cuts as in ||. Using the 
standard cuts, as derived in this paper, i.e. s min > 200 GeV and p{ ets > 7 GeV, the NLO 
and LO cross sections of MEPJET and PEPSI are compared in Table 3. For the higher order 
calculations, the largest difference is observed for the unpolarised cross sections (PEPSI/PS 
versus MEPJET/NLO) while the polarised cross sections are very similar. 
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Table 2: Asymmetries for MEPJET (LO) and PEPSI (no parton showers) for different values 
of the s cut. 
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Table 3: Asymmetries for MEPJET (LO/NLO) and PEPSI (with/without parton showers). 
The errors on the PEPSI results are of the order of 10%. 

3 Correlation of 'true' and 'visible' variables 



A problem of MEPJET is that the 'true' x g , i.e. the x value of the gluon probed in the proton, 
is not known anymore at running time, but only the x value reconstructed from the mass of the 
two jets Sij : Xj ets . In contrast to the LO case, in NLO the information of the di-jets alone is not 
sufficient to reconstruct directly x g event by event, but needs to be 'unfolded'. To be sensitive 
to AG by measuring di-jet asymmetries only there has to be a good correlation between x g 
and Xj e ts- This correlation has been checked using the program DISENT fl5|| . A correlation 



matrix has been produced, which in a second step can be used in MEPJET to reconstruct 
x g from Xj e t s . DISENT contains only unpolarised di-jet cross sections therefore we could not 
perform the whole study with this program. For the simulation we used the same conditions 
as in MEPJET, concerning jet algorithm, parton distributions and cuts. 

Fig. ^ (right) shows the correlation between x g on the x axis and Xj e ts on y-axis. The 
correlation looks promising. This correlation matrix was then applied in the MEPJET program: 
for each 'event' the x g was determined from the Xj ets randomly according to the probabilities 
of this matrix. Fig. ^| (left) shows the polarised cross sections for gluon induced processes as 
a function of Xj e ts and the so-determined x g . We see a shift to higher x as expected, but the 
corrections are not very large. 

Figure |4] shows that also in NLO the polarised cross-section for di-jet production is domi- 
nated by gluon initiated processes. The corresponding asymmetries for events due to quark and 
gluon initiated processes, which can not be distinguished on an event by event basis, are shown 
in Fig. |5] for a luminosity of 200 pb~ x . The figure also shows the asymmetries when calculated 
versus the 'measured' x (xj ets ) or the 'true' x (x true = x g ) when using the unfolding matrix 
from DISENT. It demonstrates that the effects are small due to the locality of the Xj e t s — %g 
correlation. 
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Figure 3: Left: the polarised cross section versus x (MEPJET) in LO and NLO. For the latter 
it is shown both as function of the x g and Xj ets using DISENT for the correction (see text). 
Right: the correlation between x g and Xj ets using the DISENT program. 
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Figure 4: The polarised cross section versus x in LO (left) and NLO (right) for quark and gluon 
induced processes. 



A comparison of the asymmetries expected in LO and in NLO, for the latter shown both 
as a function of x g and Xj ets , is given in Fig. [| The reduction of the asymmetry from LO to 
NLO is clearly visible. 
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Figure 5: Top: the asymmetries for gluon and quark induced processes separately versus Xj ets 
and x g , labeled as x trU e in the figure; Bottom: a direct comparison of the asymmetries versus 
Xj e t s and x g for the gluon and quark induced processes separately. 



4 Sensitivity to xAG(x) 

In this section we will quantify the sensitivity to the shape of xAG(x), following the method 
used in [[J. In a real measurement one could obtain xAG(x) from the measured asymmetry by 
an unfolding method, where the background would be subtracted statistically and correlations 
between bins are fully taken into account. The HI experiment has already shown a NLO 
extraction of the gluon by using combined information of the total inclusive and di-jet cross 
section |16| . If correlations between bins are small one can use a simpler method performing a 



bin-by-bin correction. For our study we consider the latter method to be sufficient. 

Taking the NLO GS-A gluon distribution as a reference, we calculated the statistical errors 
of xAG(x) in the range 0.005 < x < 0.4 where a significant measurement can be made. Note 
that this range is shifted slightly to higher x values compared to the LO study since x g > Xj ets . 
Also shown is the expectation for the NLO GS-C distribution. The results are shown for 
two values of the integrated luminosity and taking the polarisation for both beams to be 0.7. 
Clearly, even for a luminosity of only 200 pb~ x already a clear difference between the two gluon 
scenarios is expected. 
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Figure 6: Comparison of the asymmetries in LO and NLO for the 'true' x g and 'measured' 
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Figure 7: The statistical precision of a measurement of xAG(x) from di-jets in NLO, shown on 
top of the GS-A parton density curve, for two values of integrated luminosity The expected 
value for the beam polarisation is taken into account. GS-C is shown for reference. 



5 HERA-TESLA 



Finally the di-jet asymmetry was calculated for a possible future high energy ep collider, con- 
sisting on the one hand of the HERA proton ring, and on the other hand of a e + e~ linear collider 
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Figure 8: Asymmetries measured using a 800 GeV e~ beam of TESLA on a 820 GeV p beam 
of HERA, for two selected Q 2 regions (top/bottom). On the left hand side the low x region is 
expanded and the newly reachable low-x domain is shown by the hatched region of the plot. 



(LC). DESY proceeds towards a proposal for such a linear collider, which would have a centre 
of mass system (CMS) energy of 0.5 — 1. TeV. It is planned to include the possibility to perform 
ep collisions, by constructing the LC tangential to HERA, allowing for an interaction region in 



the HERA West hall. The kinematics and beam dynamics have been discussed in [jTT|, |I8||. The 



polarisation of the electron beam would be sufficiently large (about 80%). If also the proton 
beam is polarised, polarised ep scattering can be studied at a CMS of about 1 TeV, allowing 
to study the polarised parton distributions at an order of magnitude lower in x compared to 
HERA. In fl9| the gain for gi is discussed. Here we show the asymmetries (in LO) for the di- 
jets, using the same jet selection criteria as used for the HERA study, for collisions of 820 GeV 
protons on 800 GeV electrons, possibly the maximum which can be expected. The error bars 
correspond to 1000 pb _1 , but the polarisation of the beams is assumed to be 100%. Events and 
jets are selected within the pseudorapidity range from —4 to 3.5 for the jets and from —7 to 3.5 
for the scattered electron, Q 2 > 1 GeV 2 , pl ets > 5 GeV, and s > 100 GeV 2 . The asymmetries 
are shown for two upper limits on Q 2 . In the figures on the right the low-x region is shown 
explicitly, and the gain in x-range with respect to nominal HERA is given by the shaded area. 
The measurement reflects the decreasing asymmetries with decreasing x. The asymmetries at 
very low x become very small. The lowest values of x where a significant measurement of the 
di-jet asymmetry can be made with TESLA-HERA will be about x = 0.0005. However a large 
statistics sample 0(1) fb -1 and an excellent control of systematic errors will be needed. 
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6 Conclusion 



The direct measurement of AG(x) via di-jet production has been studied at NLO, using the 
MEP JET program. The asymmetries are reduced with respect to the LO case, but a sufficiently 
large sensitivity to the polarised gluon distribution can be obtained in the region 0.005 < x g < 
0.4 for luminosities larger than 200 pb _1 . This measurement can be extended by roughly an 
order of magnitude to lower x with future polarised ep collisions using TESLA and HERA. 
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